Multidimensional variation for quasi-Monte Carlo
نویسنده
چکیده
This paper collects together some properties of multidimensional definitions of the total variation of a real valued function. The subject has been studied for a long time. Many of the results presented here date back at least to the early 1900s. The main reason for revisiting this topic is that there has been much recent work in theory and applications of Quasi-Monte Carlo (QMC) sampling. For an account of quasi-Monte Carlo integration see Fang and Wang (1994) and Niederreiter (1992). QMC is especially competitive for multidimensional integrands with bounded variation in the sense of Hardy and Krause (BVHK). For such integrands, over d dimensional domains, one sees QMC errors that are O(n−1(log n)) when using n function evaluations. When d = 1, competing methods are usually preferred to QMC. For even modestly large d, Monte Carlo and quasi-Monte Carlo sampling become the methods of choice. When the integrand is in BVHK, then QMC has superior asymptotic behavior, compared to Monte Carlo sampling. Therefore we may like to know when a specific function is in BVHK. Recent introductory text books on real analysis typically cover the notion of total variation for functions of a single real variable. Few of them say much about multidimensional variation. The not very recent book, Hobson (1927, Chapter 5), does include some discussion of variation beyond the one dimensional case. Discussions of multidimensional variation usually require ungainly expressions that grow in complexity with the dimension d. For this reason, many authors work out details for d = 2 and report that the same results hold for all d. Yet some results that hold for d = 2 do not hold for d > 2. For example an indicator function in two dimensions must either have positive variation in Vitali’s sense, or must have at least one input variable on which it does not truly depend. The same is not true for d ≥ 3. Similarly, if f(x) and g(x) are linear functions on the d dimensional cube, then min(f(x), g(x)) is BVHK when d = 2 but is not necessarily so when d > 2.
منابع مشابه
Multidimensional quasi-Monte Carlo Malliavin Greeks
The aim of this paper is extensively investigate the performance of the estimators for the Greeks of multidimensional complex path-dependent options obtained by the aid of Malliavin Calculus. The study analyses both the computation effort and the variance reduction in the Quasi-Monte Carlo simulation framework. For this purpose, we adopt the approach employed by Montero and Kohatsu-Higa to the ...
متن کاملEfficient Bidirectional Path Tracing by Randomized Quasi-Monte Carlo Integration
As opposed to Monte Carlo integration the quasi-Monte Carlo method does not allow for an error estimate from the samples used for the integral approximation and the deterministic error bound is not accessible in the setting of computer graphics, since usually the integrands are of unbounded variation. We investigate the application of randomized quasi-Monte Carlo integration to bidirectional pa...
متن کاملACTA UNIVERSITATIS APULENSIS No 14/2007 A MIXED MONTE CARLO AND QUASI-MONTE CARLO SEQUENCE FOR MULTIDIMENSIONAL INTEGRAL ESTIMATION
In this paper, we propose a method for estimating an sdimensional integral I. We define a new hybrid sequence that we call the H-mixed sequence. We obtain a probabilistic bound for the H-discrepancy of this sequence. We define a new estimator for a multidimensional integral using the H-mixed sequence. We prove a central limit theorem for this estimator. We show that by using our estimator, we o...
متن کاملAnalysis of the Quasi-Monte Carlo Integration of the Rendering Equation
Quasi-Monte Carlo integration is said to be better than Monte-Carlo integration since its error bound can be in the order of O(N (1 )) instead of the O(N 0:5) probabilistic bound of classical Monte-Carlo integration if the integrand has finite variation. However, since in computer graphics the integrand of the rendering equation is usually discontinuous and thus has infinite variation, the supe...
متن کاملQuasi-monte Carlo Integration over a Simplex and the Entire Space
Monte Carlo integration is a widely used method to approximate high dimensional integrals. However, the randomness of this method causes the convergence to be very slow. Quasi-Monte Carlo integration uses low discrepancy sequences instead of pseudorandom sequences. The points from these sequences are more uniformly distributed. This causes the convergence to be much faster. Most research in qua...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2004